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Heat capacity measurements on the prototypical 
three-dimensional random-field Ising model compound 
Feo.58Zno.42F2 have been performed in both of the field- 
cooling (FC) and zero-field-cooling (ZFC) procedures. There 
is no evidence of hysteresis between the FC and ZFC protocols 
within the experimental accuracy. The absence of hysteresis 
is interpreted based on a magnetic field (H) — temperature (T) 
phase diagram in which a metastability region exists between 
the Neel temperature line Tn(H) and He (T), the latter be- 
ing the lowest critical field line for the local spin flips. 
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The random-field Ising model (RFIM) has been a sub- 
ject of extensive studies over the last two decades. It 
has been predicted theoretically jlj that the RFIM can 
be realized in diluted uniaxial antiferromagnets under an 
applied magnetic field. Since then, a number of exper- 
imental studies has been done on this topic. Among 
them, the three-dimensional (3D) RFIM has attracted 
much attention. The 3D RFIM shows distinctly different 
behavior depending on whether the measurements are 
done with cooling the sample in a magnetic field (field- 
cooling, FC) or done in a magnetic field with increasing 
temperature after having cooled it in zero field (zero-field- 
cooling, ZFC) below the transition temperature. Early 
neutron scattering experiments on the 3D RFIM showed 
that a long-range magnetic ordering was established in 
the ZFC case, whereas the system was in a disordered 
domain state in the FC case 

It turns out that the phase transition of the 3D RFIM 
is more complex. Jaccarino et al. jl have proposed 
a temperature (T) versus magnetic field (H) phase di- 
agram of the 3D RFIM based on experimental facts. 
They showed that a metastability region exists in the 
H — T plane bordering on the second-order phase transi- 
tion temperature line T/v(H). From a synchrotron mag- 
netic X-ray scattering study of the 3D RFIM, Hill et al. 
showed that the long-range-order (LRO) established 
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in ZFC procedure vanishes continuously with increasing 
temperature in the vicinity of T^{H). They labeled this 
as "trompe 1'oeil" critical behavior. There has been a big 
debate p|-|T^] on the nature of the phase transition of the 
3D RFIM in a recent issue of this Journal. So, the topic 
may be interesting to broad audiences. 

The phase transition in the prototypical examples of 
the 3D RFIM Fe x Zni_ x F2 has been studied by indirect 
and direct heat capacity measurements [[l3| -ftBt . From op- 
tical linear birefringence (An) measurements, Belanger 
et al. |l3f have observed a sharp peak in the temperature 
dependence of d(An)/dT of a Feo.53Zno.47F2 sample in 
finite magnetic fields (H < 2 T). Later, Dow and Be- 
langer ]l4j] have measured the heat capacity (C p ) of a 
Feo.46Zno.54F2 crystal using a variation of the classical 
heat pulse technique in H up to 1.5 T. They found an 
apparent rounding of the peak in C p for the FC measure- 
ment. This rounding of C p for the FC procedure has also 
been observed in the An study |l5| up to 1.9 T. Birgeneau 
et al. jl6| have measured C p of a Feo.5Zno.5F2 crystal us- 
ing a standard semiadiabatic dc calorimeter. They found 
no difference in C p for FC and ZFC protocols in H=1.5 
and 5.5 T. We have made a comprehensive study of C p , 
using a relaxation method, to clarify this apparently con- 
tradictory situation. In this paper, we report the results 
obtained in a single crystal of Feo.5sZno.42F2. 

The non-diluted compound FeF2 has the rutile-type 
crystal structure — PA/mnm JTtJ and it establishes 
an antiferromagnetic LRO below the Neel temperature 
Ttv=78.4 K fjqj . The spin easy axis is parallel to the c 
axis. A large single ion anisotropy makes this an excel- 
lent example of a 3D Ising system. The single crystal of 
Feo.58Zno.42F2 used in this study had been grown at the 
University of California, Santa Barbara. 

The heat capacity was measured using a MagLab HC 
microcalorimeter of Oxford Instruments, UK. This mi- 
crocalorimeter consists of a small sapphire chip on which 
a serpentine metallic heater is evaporated. Attached to 
the chip with 50 /im gold leads is a tiny temperature 
sensor. The chip is suspended by 20 /im tungsten leads 
which make electrical connections to the elements and 
also provide a weak thermal link to a calorimeter cell. 
The cell in turn is screwed onto the 3 He pot with ther- 
mally conductive grease. For a given cell temperature, 
the chip temperature is a function (only) of the power 
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dissipated in the heater element. For a measurement, the 
heater power was first increased stepwise (typically ~10 
fiW near the Neel temperature) and maintained for a pe- 
riod of 3r to ensure equilibrium was closely approached, 
where r is the relaxation time (typically ~70 sec near the 
Neel temperature). Then the power was stepped back to 
its original value. The variation of the chip tempera- 
ture with time was fitted with an exponential function 
from which r was obtained. We repeated this proce- 
dure several times to improve the quality of the data. 
With this calorimeter, FC and ZFC measurements can 
be done with the same accuracy. The single crystal was 
cut into a platelet parallel to the c axis with the dimen- 
sions about 1.5 mmx3 mmxO.l mm. The platelet was 
mounted on the sapphire chip using a small amount of 
thermally conductive grease. The heat capacity of the 
sample was obtained by subtracting the heat capacity of 
the sapphire chip from the total. The heat capacity of 
the sapphire chip is much smaller (^1/10) than that of 
the sample in the temperature range of interest. In both 
of the ZFC and FC measurements, the sample was cooled 
from T=90 K which is higher enough than Tjsr in pure 
FeF 2 . 

We measured the temperature dependence of C p of 
Feo.58Zno.42F2 in H between T and 10 T parallel to the 
c axis. Typical results under ZFC condition are shown 
in Figs.^ (a) and (b). Here, we subtracted the contribu- 
tion of the lattice from the total using C p of ZnF2 Jl8| . 
The result after subtraction corresponds to the magnetic 
part (C mag ) of the heat capacity.® In zero field (ZF) an 
asymmetric peak appears at 45.5 K. We define the peak 
temperature, T p (H), as the temperature at which C mag is 
a maximum. As is seen from Fig.[l], T p (H) decreases and 
the shape of C mag becomes symmetric and broad with 
increasing H. These observations are consistent with 
the results measured on the Feo.53Zno.47F2 (H < 2 T) 
H, Feo.45Zno.54F2 (H < 1.9 T) U) and Fe . 6 Zno. 4 F2 
(H < 8 T) [HJ samples with the An method under ZFC 
condition. 

Figure || shows the data taken at H=2 T for both of 
the FC and ZFC protocols. It is evident that there is 
no difference in the two measurements within the exper- 
imental accuracy (~0.3 K). We have observed no hys- 
teresis in other fields investigated (0.5, 1.0, and 3.0 T). 
This result is in accord with that reported by Birgeneau 
et al. |[[6| . However, the present result is in contradiction 
with the observation by Dow and Belanger [Q and by 
Ferreira et al. |lq] , in which FC heat capacity is severely 
rounded. 

These apparently contradictory results may be ex- 
plained as a result of the time scale of the measurements. 
It is widely accepted that the critical dynamics of the 3D 
RFIM is quite slow. Then, the results obtained with 
different technique such as An, adiabatic or relaxation 
method may be different. The time scale of our measure- 
ment is typically 200 sec. Moreover, we can estimate the 



equilibrium value by fitting the time evolution with an 
exponential function. 

Belanger et al. have pointed out that the concentration 
gradient of the sample makes the observation of hysteresis 
difficult @. In case of the measurement using a sample 
the large concentration gradient, it was reported that 
the peak of C mag in ZF was rounded. On the contrary, 
the peak of our heat capacity data measured in ZF is 
sharp enough to suggest that the concentration gradient 
of our sample is very small. In order to confirm this, we 
analyzed the data measured with ZF by fitting them to 
the scaling function pp[ 
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where t=T /T/v-1. The fitting parameters thus obtained 
close to those reported from the analysis of the indirect 
heat capacity data, d(An)/dT, measured in ZF using the 
sample of Feo.60Zno.40F2 |2C|| . Thus, we are convinced 
that our sample is homogeneous enough to observe the 
random field effects. 

Next, we discuss the apparent rounding of the peak at 
high fields (Fig.|| (b)). In Fig.||, we show the field depen- 
dence of the width (W p ) of the peak of C mag defined as 
the one at which the value of C mag is 95 % of its maxi- 
mum in respective H obtained from the ZFC procedure. 
As is seen from Figj^, W p is almost constant below 4 T. 
On the other hand, W p depends much on H above 4 T. 
We interpret this broadening of the peak of C mag as an 
evidence that the system is in a domain state at finite 
fields. 

King et al. |2lJ have observed in Fe 2: Zni_ :c F2 that 
the magnetization shows anomalies at He — (n/8)HE, 
where He is the exchange field and n=l, 2, ... 5. These 
anomalies have been explained as arising from local spin 
flips which occur when the external magnetic field be- 
comes equal to the local exchange fields H^ ■ The low- 
est field H^ required for the flip is about 7 T at 1.3 
K pl| . Since the exchange field is proportional to the 
thermal average of the magnetic moment, it decreases 
with increasing temperature and vanishes at TV. When 
the sample is cooled in zero field below T/v and H = Hq 
is subsequently applied, the system is in the LRO state 
provided H < He (T). At an elevated temperature, 
Hc X \t) becomes equal to Hq, where the local spin flip 
occurs. This spin flip will nucleate domains. Near the 
Neel temperature, H^ (T) with n=l, 2, ... become close, 
thereby domains with various sizes will be nucleated. It 
is not surprising that the "trompe l'oeil" critical behavior 
has been observed in the 3D RFIM at finite fields ||[l6) . 

When the sample is cooled from high temperature un- 
der a magnetic field (FC case) and T/v (if) is reached, 
a LRO will be established there. However, as is dis- 
cussed above, Hc n \T) with n=l, 2, ... are close near 
Tn(H), thence domains will be nucleated immediately 
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below Tjv(i?). These domains are stable at low temper- 
atures against thermal agitation because the anisotropy 
is strong. According to this interpretation, the transi- 
tion temperature in equilibrium is the Neel temperature 
Tn(H) and there is a metastability region below Tjv(iJ) 

bounded by the H^\t) line. This conclusion is con- 
sistent with the result of the computer calculation ]22j | 
made on a diluted 3D Ising antiferromagnet in magnetic 
fields. 

Because the system is in the domain state at this 
metastability region for both of the ZFC and FC con- 
ditions, the heat capacities coming from thermal fluctu- 
ation in these states are expected to be similar. This ex- 
plains the absence of hysteresis in the ZFC and FC heat 
capacities as observed in the present experiment. Based 
on this interpretation, then, it is possible to observe 
the crossover from the random exchange Ising model to 
RFIM with increasing H along the Tjv(-ff) line. 

In the remaining part of this paper, we discuss the shift 
in the transition temperature with H. Here, we assign 
the peak temperature T p (H) as the transition tempera- 
ture. The true transition temperature Tn(H) may be a 
bit larger than T p (H) for a given H . However, as is seen 
from Fig.||, the difference between T p {H) and Tpj(H) lies 
within the experimental error below 4 T. 

In Figg, AT p (H)(= T p {H = 0) - T p {H)) is plotted as 
a function of H on a log-log scale. Data were fitted to 
Eq.(2) below over the field range 0.5 T <H< 4.0 T from 
the reason discussed above, 
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where a and b are constants and <j) is the crossover expo- 
nent. The term bH 2 represents a small mean-field shift. 
From the fitting, we got <fi — 1.39 ± 0.03. This value is 
close to the ones observed experimentally 13 ll]|2j| and 
the theoretical ones (0=1.25 0-1.4 |J]). 

In conclusion, we have made a comprehensive study of 
C p , using a relaxation method, on the prototypical 3D 
RFIM compound Feo.5sZno.42F2 in both of the FC and 
ZFC procedures. We have observed no hysteresis in the 
two measurements within the experimental accuracy. We 
interpret this result based on an H — T phase diagram in 
which a metastability region exists between Tpj(H) and 
He (T), the latter being the lowest critical field for the 
local spin flip. This phase diagram seems to settle some 
of the conflicting experimental results reported so far. 
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FIG. 1. Temperature dependence of the magnetic heat ca- 
pacity in a Feo.58Zno.42F2 single crystal obtained under ZFC 
condition. The external magnetic field is applied parallel to 
the c axis. 



FIG. 2. Temperature dependence of the magnetic heat ca- 
pacity in a Feo.58Zno.42F2 single crystal taken at 2 T for FC 
and ZFC protocols. 
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FIG. 3. Field dependence of the width of the peak in the 
magnetic heat capacity of a Feo.5sZno.42F2 single crystal. 

FIG. 4. Shift of the peak temperature in the magnetic heat 
capacity of a Feo.58Zno.42F2 single crystal with magnetic field. 
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